Observation of diffractive orbits in the spectrum of excited NO in a magnetic field 



A. Matzkin 

Laboratoire de Spectrometrie physique ( CNRS Unite 
Universite Joseph-Fourier Grenoble-I, BP 87, 38402 Saint-Martin, France 



o 
o 

(N 
X 



in 



M. Raoult 

Laboratoire Aime Cotton (CNRS Unite 3321), Universite de Paris-Sud, 91405 Orsay, France 

D. Gauyacq 

Laboratoire de Photophysique moleculaire ( CNRS Unite 3361 ), Universite de Paris-Sud, 91405 Orsay, France 

We investigate the experimental spectrum of excited NO molecules in the diamagnetic regime and 
develop a quantitative semiclassical framework to account for the results. We show the dynamics 
can be interpreted in terms of classical orbits provided that in addition to the geometric orbits, 
diffractive effects are appropriately taken into account. We also show how individual orbits can 
be extracted from the experimental signal and use this procedure to reveal the first experimental 
manifestation of inelastic diffractive orbits. 
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Periodic orbit (PO) theory Q has been a very success- 
ful quantitative and qualitative tool in the analysis of 
systems displaying strong quantum fluctuations. In par- 
ticular, the complex photoabsorption spectrum of excited 
hydrogen atoms in a static magnetic field was interpreted 
within closed orbit theory (COT) in terms of periodic or- 
bits closed at the nucleus [2|. However standard PO the- 
ory was shown to be insufficient for systems containing a 
scatterer smaller or comparable to one de Broglie wave- 
length, because the PO's may hit the scatterer, thereby 
producing diffractive orbits In certain mesoscopic 
devices diffractive orbits crucially influence the spectrum 
whereas in non-hydrogenic Rydberg atoms in static 
fields the atomic core acts as a scatterer producing com- 
bination orbits 5]. This was experimentally observed in 
the photoabsorption spectrum of helium, which departs 
slightly from the hydrogenic case Q . In the case of Ryd- 
berg molecules, core effects are expected to be more im- 
portant: the reason is that the diffractive scatterer has 
an internal structure (the quantum states of the molec- 
ular core). Indeed, model calculations recently predicted 
that Rydberg molecules in external fields would display a 
novel type of diffractive orbit produced by inelastic scat- 
tering of the PO's on the core 0. 

In this work, we address this problem by investigating 
the experimental spectrum of excited NO molecules in 
a magnetic field. We show the apparently complex spec- 
trum can be well understood in terms of closed orbits pro- 
vided molecular core effects are appropriately taken into 
account; in particular, we shall see COT calculations cor- 
rectly reproduce the Fourier transform of the experimen- 
tal spectrum, in a regime in which quantum calculations 
have yet to be undertaken. We will also demonstrate how 
individual orbits may be extracted from the experimental 
results by using a modified Fourier tranform, and use this 
technique to show the first experimental manifestation of 



an inelastic diffractive orbit in a physical system. 

Most previous studies of atoms in a magnetic field have 
relied on scaled-energy spectroscopy, in which the field 
B is varied with the Rydberg electron energy E so as 
to keep the scaled energy e = E(B/B )- 2 / 3 constant, 
with Bq = 2.35 x 10 5 T. Then the classical motion is 
unchanged at each point of the spectrum and each orbit 
contributes as a peak in the Fourier transformed (FT) 
spectrum. However, such a scaling law does not hold, al- 
beit approximately, for molecules. We have therefore ob- 
served the Rydberg states of NO at a constant magnetic 
field of 0.93 T, varying the excitation energy of the probe 
laser. The experiment, carried out in a "magnetic bottle" 
time of flight spectrometer, was described elsewhere ||, 
but in short: a 2-photon highly selective transition with 
a pump laser brings the molecule from the ground state 
X 2 H 3 / 2 (in the vibrational level v" = 0) into the single 
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1: Experimental spectrum for the transition from the 
A 2 T, + , N = 3, Mn = —3, 3 states of NO in a static mag- 
netic field of strength B — 0.93 T. The lower scale gives the 
energy relative to the ground state of NO. The origin of the 
upper scale corresponds to the field-free ionization threshold 
for levels converging to the N + = 1 rotational state of NO + . 
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FIG. 2: (a) Fourier-transform of the experimental spectrum. 
The first peak (marked with a star) is a residual oscillation 
from the background due to the redistribution of the oscil- 
lator strength into molecular Zeeman structures at low ener- 
gies 9] . (b) Closed-orbit semiclassical calculations with (solid 
line) and without (dotted line) the diffractive contributions. 
An independent calculation was performed for each value of 
Mjv; each result was then multiplied by the relative weight 
of the Mjv component before taking the Fourier transform. 
Note that the wide peaks are not correlated with single closed 
orbits but arise from a coherent superposition of orbits at dif- 
ferent energies. The shortest contributing orbit is an orbit 
perpendicular to the field with a period of 4.22 ps. 



rotational Zeeman sublevel N = 3 Ms — 1/2 of the first 
electronic excited state A 2 £ + , v' = 1; N is the total an- 
gular momentum exclusive of spin and the constant spin 
number Ms will henceforth be disregarded. The sub- 
levels Mm = —3 to 3 (where Mn is the projection of 
N on the field axis) were not resolved but their relative 
population is known. In a second step, the probe laser, 
linearly polarized parallel to the field axis, excites the 
molecules to high-lying Rydberg states having also v = 1 
(so we can neglect vibrational couplings). 

These excited states are further ionized via vibrational 
autoionization into the ground state of the NO + core. 
The resulting spectrum, which lies higher in energy than 
those previously published in [8| (in which nonperturba- 
tive diamagnetic effects could be neglected), is shown in 
Fig. 1. The FT of the experimental spectrum is shown 
in Fig 2(a), whereas Fig. 2(b) displays the semiclassi- 
cal calculations which will be detailed below. Unlike a 
spectrum taken at constant e, we do not expect in the 
present case to correlate each peak with a single orbit, 
because the period of an orbit changes with the energy. 
Before turning to the extraction of individual orbits, we 
detail the mechanism giving rise to the peaks in the FT 
spectrum. 

Since in the A 2 T, + state the orbital momentum of the 
Rydberg electron is a mixture of both / = and 1 = 2, 
the probe laser excites to final states having I = 1,3 and 
N = 2,3,4. Photo-excitation takes place in the molec- 
ular frame (in which the Rydberg electron is coupled to 



the molecular core), but following excitation the electron 
uncouples from the core and N is no longer conserved, 
although its projection Mpj is The molecule is then 
described as an electron orbiting around the molecular 
core NO + in well-defined rotational states N + = 1,3 or 
5. For each core configuration, the electron's wavefunc- 
tion is propagated semiclassically along classical trajec- 
tories (see Fig. 3). The energy of the Rydberg electron 
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FIG. 3: Schematic diagram of the excitation dynamics of the 
Rydberg states of NO in a magnetic field. (1) A pump laser 
excites the molecules from the ground state X 2 Yl 3 ^ 2 to the 
intermediate state j4 2 E + , N = 3, via a 2-photon transition. 
A high-resolution probe laser then excites the molecules from 
the intermediate state to Rydberg levels with orbital momen- 
tum I = 1 and / = 3. (2) Following this excitation the Ryd- 
berg electron uncouples from the molecular core which is left 
freely rotating with rotation quantum numbers N + = 1, 3 
or 5. (3) As the Rydberg electron leaves the core region its 
dynamics is described using classical trajectories; the molec- 
ular energy is partitioned between the Rydberg electron and 
the core: as the core rotation increases (pictured by a larger 
core), the energy of the Rydberg electron decreases (sym- 
bolized by smaller orbits). The core region is less than 5 a.u. 
large, whereas even the shortest orbit associated with a highly 
rotating core in state N + = 5 extends beyond 1500 a.u. (4) 
Some trajectories return to the core, which following the ter- 
minology of diffractive PO theory we call the geometric closed 
orbits; they are the same orbits that exist in the diamagnetic 
hydrogen atom (i.e. in the absence of the core). (5) As the 
Rydberg electron enters the core region, its dynamics recou- 
ples with the core. The waves carried by the geometric CO's 
overlap with the initially-dipole excited waves, producing fluc- 
tuations in the oscillator strength. (6) The electron scatters 
off the core, uncouples, and follows again classical trajecto- 
ries; during the scattering process, the electron may have ex- 
changed energy with the core. Eventually, some trajectories 
return once more to the core; these diffractive orbits, com- 
posed of 2 geometric CO's connected by elastic or inelastic 
core-scattering, produce further fluctuations in the oscillator 
strength. 
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-1.1 — > -0.4 
0.133 
0.004 


-1.5 — > -0.7 
0.304 
0.264 


-2.1 — > -1.4 
0.225 
0.003 



TABLE I: For each state of the core, we give the range of the 
scaled energy e of the Rydberg electron as the probe laser is 
scanned. When N + = 5, the very low values of e correspond 
to the near integrable regime. For N + = 3, the dynamics 
of the Rydberg electron is mostly regular. For N + = 1, the 
electron enters the chaotic regime: at e = —0.4, about 50 % 
of the classical phase space is chaotic |l2^1 . We also give the 
modulus of Q° and Q D for N +l = 3 and for M N = -1; these 
values give a rough idea of the PO-independent part of the 
fluctuations strength associated with the different geometric 
and diffractive contributions. 
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=1 +B r [N+(N+ + l)-2]+^B/B , (1) 



where B r is the rotational constant of NO + . Disregarding 
the linear Zeeman shift, we have three main dynamical 
regimes for the Rydberg electron, each associated with 
a different rotational core state. Table I gives, for each 
value of N + the scaled energy ranges corresponding to 
the experimental spectrum. 

At each energy E N +, some of the trajectories are 
turned back by the field and return to the core; those 
trajectories are the geometric closed orbits. Their contri- 
bution to the oscillator strength f(E) is given by |9j 
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Ah and Sk are 3-dimensional quantities respectively 
linked to the 2-dimensional classical amplitude and 
action Sk of the fcth closed orbit ^lj (<->fe includes the 
Maslov index fj,k(Ej^+)); those quantities depend on the 
energy of the Rydberg electron, i.e. on the core state and 
on the excitation energy. <2^+ is a quantity depending 
on the dipole transition factor from the A 2 T, + to the ex- 
cited states, and on the transformation elements from the 
coupled to the uncoupled frame; some values are given in 
Table I. 

When the semiclassical wavefunction carried by the ge- 
ometric closed orbits returns to the core, it diffracts on 
the scatterer (Fig. 3). This diffraction process, produces 
newly outgoing trajectories q, which eventually return to 
the core; the resulting diffractive contribution to f(E) is 
given by 

E Q-N+ ,N+' ^ ( E N+ )-A q (E N +, ) 
N+,N+' k,q 

x expi [Sk{E N +) +S q (E N +,)\ . (3) 



Qn+ jv+' acts as a global diffraction coefficient; unlike 
diffraction in quantum billiards 0,0], Q%+ N +, does not 
depend here on the local geometry of the scatterer but 
solely on the quantum properties of the core (i.e. on 
the core-induced phase-shifts, known as molecular quan- 
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FIG. 4: Extraction of individual orbits from the experimen- 
tal spectrum using the chirped FT procedure. In the present 
illustration, we show the results for 4 orbits which exist at 
an excitation energy Eq — —26.91 cm - below the held free 
N + — 1 ionization threshold, (a) Extraction of the second re- 
turn of the perpendicular orbit (R2 in the nomenclature given 
e.g. in [hUp. associated with the core in state N + — 5; the 
corresponding scaled energy of this orbit is e = —1.49. Broken 
curve: standard FT of the experimental spectrum, shown in 
Fig. 1. Thick curve: chirped FT of the experimental spec- 
trum obtained by calculating the function £(i; E) appearing 
in Eq. @ for this orbit. The arrow indicates the calculated 
classical period of R2, which agrees very well with the posi- 
tion of the peak appearing in the CFT. (b) CFT extraction 
of the orbit Vi, parallel to the field axis, associated with the 
core in state N + — 3, at e = —0.85. (c) Extraction of the Rl 
orbit associated with the core in state N + = 1; orbits asso- 
ciated with N + — 1 are generically more difficult to extract 
from the experimental signal because the oscillator strength 
tends to be dominated by orbits associated with N + — 5 
and 3 (see Table I). However, at this energy, corresponding 
to e = —0.48, the classical amplitude A R i is very strong (be- 
cause e is close to the bifurcation energy at which this orbit is 
created [TBI'), thereby allowing the signal to pick up its con- 
tribution, (d) By applying the CFT procedure to the action 
expansion Sv ± (E N + =3 ) + S R i(E N + =1 ), we can also extract 
the diffractive orbit formed by the combination of the orbits 
extracted in (b) and (c); the period of the diffractive orbit 
is given by summing the periods of the two geometric CO's. 
The shape of the orbits is given in the upper right-hand side 
of each plot (B is in the vertical direction). 
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turn defects; the value of the 6 independent quantum de- 
fects for NO is given in 8]). The matching between the 
semiclassical waves and the core wavefunction takes place 
in the stationary phase approximation. The oscillator 
strength contains both elastic (N + = N + ') and inelastic 
(N + ^ N +l ) diffractive terms. In the one core-scatter 
approximation, the total contribution to the fluctuations 
in the oscillator strength is obtained by summing the 
imaginary parts of Eqs. (0 and ©. 

The experimental or semiclassically calculated FT 
spectra are respectively obtained by Fourier 

transforming the photoelectron signal or the oscillator 
strength for which we shall use the same notation f(E), 
hence 

T(t; = [ ^ f{E) exp E)) dE (4) 

where £(t,E) = Et, and we set the energy variable rela- 
tive to the N + — 1 ionization threshold as E = E N + =1 . 
The semiclassical calculations were undertaken by apply- 
ing Eqs. J2J and Since at each point E of the spec- 
trum the orbits and therefore their amplitude and action 
change, we have determined the classical quantities A k , 
Sk, /ifc on an energy grid for the primitive closed orbits 
and their repetitions, keeping track of bifurcations as the 
energy increases. We then interpolated those functions 
so as to obtain classical quantities varying with the en- 
ergy. This procedure, which is tractable for short orbits 
provided the scaled energy is not too high, yields a good 
agreement with the experimental FT spectrum, as seen 
on Fig. 2(b). The semiclassical calculations allow to sin- 
gle out contributions from the ensemble of geometric or 
diffractive orbits associated with given core states. Full 
details will be given elsewhere @ . 

Eq. (0} is also at the basis of the extraction of indi- 
vidual orbits from the experimental spectrum provided 
£ (t, E) is appropriately modified so as to counterbalance 
the varying values of the classical quantites. Indeed in 
the neighborhood of some energy Eq a classical orbit k 
contributes locally as a sinusoidal oscillation of period 
Tk(Eo). But as the energy is varied, the classical motion 
changes, turning the sinusoidal oscillation into a chirped 
modulation. It has been shown very recently that a 
chirped Fourier Transform (CFT) converts a chirp mod- 
ulation into a sine oscillation, thus giving a sharp peak in 
the FT spectrum. A CFT is defined by setting £(t, E) = 
t(E - E ) + b{E - E Q ) 2 + c(E - E ) 3 + i\na(E - E ) 
in Eq. b and c are numerical coefficients and 

a is a function, which depend on the orbit to be ex- 
tracted: to extract an orbit k at some energy Eq, we 
Taylor expand the action Sk(E). To first order, we have 
dES(Eo) = Tk(E ), i.e. the period of the orbit. There- 
fore by setting b = -d%S k (E Q )/2 and c = -d%S k (E Q )/6 



in the CFT exponent £, we can expect to effectively lin- 
earize the signal around Eq so that the CFT yields a 
peak at t = Tk(Eo). The peak can be enhanced by us- 
ing a function a to compensate for the signal amplitude 
variation; we have used the inverse of our interpolated 
functions Ak(E). As discussed in great details in [l3| . 
spurious peaks inevitably appear in the CFT spectrum. 
We proceed in the following manner: using purely clas- 
sical calculations, we compute on the one hand the CFT 
of the experimental signal with the relevant parameters 
a, 6, and a, and on the other hand we determine the 
period Tk(Eo); the consistency of the procedure is en- 
sured by the appearance of a peak at the computed value 
t = Tk(Eo) when the orbit is succesfully extracted from 
the signal. An illustration of the extraction procedure is 
given in Fig. 4. For a fixed energy Eq we extract one CO 
associated with each core state (panels (a)-(c)). Panel 
(d) shows the extraction of the diffractive CO formed by 
the combination of the orbits extracted in (b) and (c) . 

In summary, we have investigated the experimental 
spectrum of NO in the diamagnetic regime and inter- 
preted the results in terms of geometric and diffractive 
classical orbits, and we have shown how individual orbits 
can be extracted from the experimental signal. 
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